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A Bloch Point singularity can form a metastable state in a magnetic nanosphere. We classify 
possible types of Bloch points and derive analytically the shape of magnetization distribution of 
different Bloch points. We show that external gradient field can stabilize the Bloch point: the shape 
of the Bloch point becomes radial-dependent one. We compute the magnetization structure of the 
nanosphere, which is in a good agrement with performed spin-lattice simulations. 
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INTRODUCTION 

Topological singularities are widely recognized as key 
to understanding the behavior of wide variety condensed 
matter systems. Linear topological singularities such as 
dislocations, disclinations, and vortices, play a crucial 
role in low-dimensional phase transitions,^ crystalline 
ordering on curved surfaces^, rotating trapped Bosc- 
Einstein condesateiP etc. Recent advances in micro- 
structuring technology have made it possible to fabri- 
cate various nanoparticles with well-prescribed geome- 
try. Much recent research in this field has focused on 
the statics and dynamics of topological singularities in 
nanoscale confined systems: essentially inhomogeneous 
states can be realized in magnetic nanoparticles-^ and 
ferroelectric nanoparticled^. As a result of the compe- 
tition between exchange and magnetic dipole-dipole in- 
teractions the ground state of magnetic disks with sizes 
larger than some tens of nanometers is a flux-closure vor- 
tex state. 

Besides linear singularities there exist also so-called 
point singularities such as monopoles, Bloch points, boo- 
jums. For example, hedgehog (monopole) singularities 
play a crucial role in the behavior of matter near quantum 
phase transitions that are seen in a variety of experimen- 
tally relevant two-dimensional antifcrromagnetsP boo- 
jums are relevant in superfluid He— 3,^ Bloch points along 
with Bloch lines are p rinci ple in understanding of mag- 
netic bubble dynamics 

The concept of point singularities was introduced in 
magnetism by FeldtkellerEU, who considered different 
magnetization distributions inside the singularity and 
proposed first estimations of the Bloch point shape. 
Later Doring^ studies how magnetostatic energy gov- 
erns the Bloch point structure by selecting the rotation 
angle inside the Bloch point. Bloch point singularities 
were directly observed in yttrium iron garnet crystals.^ 
During the last decade Bloch points were also studied 
by micromagnetic simulations in na nowiresj^ in bub- 
ble materials,^ in disks-shape d 16 * 17 -* and astroid-shaped 
nanodots^. The ultrafast switching of the vortex core 
magnetization open doors to consider the vortex state 
nanoparticles as promising candidates for magnetic ele- 



ments of storage devices. There are different scenarios 
of the switching process: (i) The symmetric or so-called 
punch-through core reversal takes place under the ac- 
tion of DC m agnetic fie ld applied perpendicularly to the 
magnet planeP^ 16 * 19 * 20 * This reversal process as a rule is 
mediated by creation of two Bloch points.^ However sin- 
gle Bloch point scenario was also mentioned in Thiaville 
et alJ^. (ii) The switching under the action of differ- 
ent in-plane AC magnetic fields or by a spin polarized 
currents pIHUD [ s accompanied by the temporary creation 
and annihilation of the vortex-antivortex pair. The lat- 
ter is accompanied by Bloch point creation^. 

The purpose of the current work is to study the mag- 
netization structure of the Bloch point of the spherical 
nanosized particle. As opposed to bubble films, where 
the static Bloch poin t results from the transition be- 
tween Bloch linespESl and vortex nanodots, where the 
Bloch point dynamically appears during the vortex core 
switching processp^i' the Bloch point in the nanosphere 
is an example of "pure" singularity without surrounding. 
Such a singularity is in some respect the only stable sin- 
gularity in ferromagnetic We consider different types of 
Bloch point and classify them in terms of vortex parame- 
ters. The conventional magnetization distribution in the 
Bloch point is generalized for the radial-dependent one. 
Such radial distribution becomes important for the Bloch 
point nanosphere under the action of nonhomogeneous 
magnetic field. We show that radial gradient field can 
stabilize the Bloch point and compute the magnetization 
structure, which is in a good agrement with performed 
spin-lattice simulations. 

The paper is organized as follows. In Sec. |T] we de- 
scribe the model and present the classification of differ- 
ent Bloch point types (Sec. I A I. The energetic analysis 
and the Bloch structure is analyzed in Sec. |IB| In order 
to stabilize the Bloch point inside the nanosphere, we 
consider the the influence of external gradient field on 
the magnetization structure. The Bloch point solution 
becomes radially dependent: we calculate the magnetiza- 



tion structure analytically in Sec. [II] In Sec. Ill we study 
the Bloch point structure numerically, in particular, the 
problem of stability. We discuss our results in Sec. |IV| In 
Appendix |A"| we analyze the Bloch point structure under 



2 



the influence of weak fields using the linearized equations. 



THE MODEL AND THE BLOCH POINT 
SOLUTIONS 



The equilibrium magnetization configuration is deter- 
mined by minimization of the energy functional ([I]), 
which leads to the following set of equations: 



eV 2 m = Vip, VV = V 



(4) 



Let us consider the classical isotropic ferromagnetic 
sphere of the radius R. The continuum dynam- 
ics of the magnetization can be described in terms 
of the magnetization unit vector m = M/M$ — 
(sin 8 cos $, sin O sin cos O) , where 9 and <& are, in 
general, functions of the coordinates and the time, and 
Ms is the saturation magnetization. The total energy E 
of such a sphere, normalized by AttM^V with V = §7T-R 3 
reads: 



S = £ ex + S l + <r 



(la) 



The first term in (la) is dimensionless exchange energy 



8tt 



dr 



(V9) z + sin-'e(V$)^ 



(lb) 



with e = £ 2 /R 2 being reduced exchange length, I — 
\J A i '4m Ma being the exchange length, A being the ex- 
change constant and r — (x,y,z)/R being the reduced 
radius-vector. The second term determines the interac- 
tion with external magnetic field H: 



3 

4-7T 



J dr (m ■ h) 



(lc) 



where h = H/AttMs is a reduced external field. We will 
discuss the influence of external field later, see Sec. [IT] 
The last term determines the reduced magnetostatic en- 
ergy: 
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dr (m ■ h r 



(Id) 



where h ms = H ms /4-kMs is a reduced magnetostatic 
field H ms . Magnetostatic field h ms satisfies the Maxwell 
magnetostatic equations^' 



V x h ms = 0, 

V • h ms = 4ttA, 



(2) 



which can be solved using magnetostatic potential, 
h ms = —"Vip. The source of the field h ms are magne- 
tostatic charges: volume charges A = —(V • m)/47r and 
surface ones a = (m ■ n)/4ir with n being the external 
normal. The magnetostatic potential inside the sample 
reads: 



ip(r) 



dr 



, A(r') 



dS'j^- (3a) 



— f dr'(m(r')-V r ,) - — — r . 
4tt J V v ' r J |r - r'| 



(3b) 



A. Classification of singularities 

Let us start the Bloch point as a particular solution of 
Q. In the exchange approach the simplest hedgehog- 
type Bloch point is characterized by the magnetization 
distribution of the form m = r/r with a singularity at 
the origin. Using spherical frame of reference for the 
radius-vector r with the polar angle i? and azimuthal 
one (p, one can describe the magnetization angles of such 
a Bloch point as follows: 9 = i) and $ = ip. The energy 
of the Bloch point in the exchange approach reads^ 



= AttAR. 



(5) 



This interaction is invariant with respect to the joint ro- 
tation of all magnetization vectors, which gives a possibil- 
ity to consid er family of solutions with different rotation 
angles .H1I12I 

We consider the following singular magnetization dis- 
tribution: 

9(tf) =ptf+n(l-p)/2, $(<p) = qtp+-f, p,q = ±l, (6) 

which describes a three-parameter Bloch point. We re- 
fer to the parameter q = ±1 as to vorticity of the Bloch 
point and p = ±1 as to its polarity using the conven- 
tional symbols for magnetic vortices. The last parameter 
7 des cribe s the azimuthal rotational angle of the Bloch 
point™ 

We refer to the micromagnetic singularity Q as to 
BPP. For example, the hedgehog-type Bloch point is a 
vortex Bloch point with positive polarity (p = 1, q = 1, 
7 = 0). The schematic of magnetization distribution in 
different types of Bloch points is presented on Fig. [TJ The 
analogy between Bloch point and vortices comes from the 
symmetric or punch-through vortex polarity switching 
process under the action of DC perpendicular magnetic 
fieldPThis reversal process as a rule is mediated by cre- 
ation of two Bloch points!^ For example, two singulari- 
ties, BP} and BP] -1 describe intermediate state between 
two vortices with opposite polarities, see Fig. [lc] It is 



also instructive to mention that a single Bloch point can 
be imagined as a composite of two vortices with oppo- 
site polarities: such a singularity can appear in 3D Eu- 
clidean space durin g the vortex polarity switching process 
in antifcrromagnets^H. All four distributions for differ- 
ent signs of p and q can be observed during symmetrical 
Bloch poi nts injection in polarity switching process of 
vortices^ and antivortices^H 

Topological properties of the Bloch point can be de- 
scribed by the topological (Pontryagin) index 



-in 



sin9(r)d9(r)d$(r) = pq. 



(7) 
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(a) p = q = 1 (b) p = — 1 q = 1 (c) Vortex (d) p = 1 g = -1 (e) p = g = — 1 (f) 

switching Antivortex 

switching 



FIG. 1: Schematic of different types of Bloch points. Magnetization distribution in azimuthal vortex Bloch points in 
sphere, see Figs, la lb and both Bloch points in axial part of cylinder-shaped sample during the vortex polarity 
switching process, see Fig. [lc] The same is for azimuthal antivortex Bloch points, see Figs. |ld[ Fig. [TeJ and both 
singularities in axial part of astroid-shaped sample during the switching, see Fig. [Tfj 



Different Bloch point distributions with equal Q are topo- and obtained 
logically equivalent: e.g., BPl} can be obtained from 
BP} by simultaneous rotation of all magnetization vec- j D 
tors by tt in vertical plane, and BP^i transforms to BP^ 1 
by rotation by it/2 in vertical plane. Note that similar 
topological notations were introduced by Malozemoff and 
Slonzewski^ for magnetic bubbles!^ 
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B. Magnetization structure of Bloch points 

The most strong exchange interaction is invariant with 
respect to the rotation angle 7. Such degeneracy is re- 
moved under account of magnetostatic interaction. It is 
worth noting that the problem of stray field influence on 
the Bloch point energetics has a long story. Feldtkeller in 
his pioneer work^ used a so-called pole avoidance princi- 
ple, see e.g. Ref.l2"9"l the magnetostatic tries to avoid any 
sort of volume or surface charge. In this way he calcu- 
lated the angle 7 from the condition that the total volume 
magnetostatic charge J X(r)dr = 0, where A(r) is the 
charge density. For the Bloch point given by Ansatz Q 
it has a form A(r) = — [psin 2 $ + cos7(cos 2 i9 + 1)] /47rr 
and leads to the rotation angle 



7f 



H) 



120° p = +1, 
60° p = -1. 



(8) 



In it interesting to note that the same value 7f also 
corresponds to absence of the total surface charge, 
/ a(r)dS = 0, where the surface charge density a{r) — 
(p cos 2 z9 + cos 7 sin 2 /4tt. 

Another approach was put forward by Doring^, who 
determined the equilibrium angle of 7 by minimizing the 
energy 

'* 1 dr(h ms ) 2 (9) 



8tt 



However one has to emphasize that the equilibrium angle 
( 10 1 minimizes only the inner part of the magnetostatic 
energy because the integration in ^ is carried over the 
sample volume V while the outer part of stray field is 
ignored. Note the similar approach was used in quite 
recent paper p^l where a magnetization contraction was 
taken into account. 

The aim of this section is to find the equilibrium rota- 
tion angle which minimizes the total magnetostatic en- 
ergy. In order to derive the magnetostatic energy of 
Bloch points ([6]), we calculate first magnetostatic poten- 
tial ( 3b I using an expansion of l/\r— r'\ over the spherical 
harmonics, 



1 

|r — r' 



00 1 



;=0m=-Z 



47T 

27+lVr 



with r< = min(r, r') and r> = max(r, r') which results 



^0=1 ( r ) — V^t H — (9 r — 8) cos 7 + nr (p — cos 7) cos 2 
3 

■0g = _ 1 (r) = p7rr(l + cos 2 1}) + irrcos{2ip + 7) sin 2 1). 

Simple calculations show that the magnetostatic energy 
of the antivortex Bloch point does not depend on 7 and 
<g'™_ 1 = 7/30 « 0.23. In contrast to this, the vortex 
Bloch point energy depends on the rotation angle 7 and 
has the form 

<Ti P (7) = ^ (7 + 4p cos 7 + 4 cos 2 7 ) . (12) 

The equilibrium value of rotation angle 70 corresponds 
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FIG. 2: (Color online.) The Bloch point energy vs 
rotation angle for BP J: analytical result |l2| (solid 

curve) and simulations (symbols). Simulations 
parameters: sphere diameter 2R — 35ao, exchange 

length £ — 3.95ao, damping parameter r\ = 0.5. 



to the minimum of the energy (12). It gives 
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105°p = +l, 
76° p = -l. 



(13) 



the question: How to stabilize the Bloch point? In this 
section we show that one way to achieve this goal is to 
apply a magnetic field which has the same symmetry as 
the hedgehog Bloch point with m = r/r, i. e. a radial 
symmetric magnetic gradient magnetic field in the form 



h = br 



(15) 



Under the action of the space dependent magnetic field 
( |15| the magnetization distribution also becomes space 
dependent. We take into account possible dependence by 
the following radial Bloch point Ansatz 

8(i?) = p-d + ?r(l - p) /2, $(r, tp) = qip + 7(r) (16) 

with a radially dependent parameter j(r) in comparison 
with Eq. ([6]). The form of this Ansatz will justified by 
numerical simulations in Sec. Mil 



Inserting Eq. ( 16 1 into Eq. ( lb ) for the exchange en- 



ergy of such magnetization distribution we get 



3e + e 



dr 



r dr. 



(17a) 



The magnetostatical potential of the Bloch point ( 16 ) 
reads 



Let us compare Bloch point energies ( 12 ) for above 
mentioned approaches: the energy of Feldtkeller 11 Bloch 
point ^? 9 ™ s 1 p (7f) = 0.1, for D6ringQ21 Bloch point one 
has # ms (7 D ) « 0.088, the result by Elfas and Verga^ 
is <^' 1 ms (7Ev) ~ 0.089. The minimal energy has a Bloch 
point with the rotation angle 70, see (13): 



C7(7o) = 



12 



0.083. 



(14) 



In order to verify our results we performed numerical 



spin-lattice simulations, see details in Sec. Ill We com- 
pare analytical dependence <CI1V _1 ( 7 ), see Eq. ( 12 ), with 



the discrete energy p4l , extracted from simulations, see 
Fig. [2j Both depende ncies are matched in maximum at 
7 = 0. Comparison can be provided by calculating en- 
ergy gain A J 5 (7) = - ^(j) for different rotation 
angles 7. According to simulation results the energy gain 
for mentioned above angles read: 

A<? (tf) ps 0.446, A<?( 7D ) w 0.460, AS" ( 7o ) ps 0.465. 

The maximum energy gain takes place for 70, which cor- 
responds to the energy minimum in a good agrement with 
our analytical result (13). 



II. THE BLOCH POINTS IN EXTERNAL FIELD 

The Bloch point does not form a ground state of a 
magnetic sphere. It corresponds to the saddle point 
(sphaleron) of the energy functional. This brings up 



<Z» = -y / ! +2<-(,s -(,'.)] <1,'- 



4tt 3 cos 2 - 1 



r' 3 [cos 7(7-') - 1] dr'. 



Here and below we consider the case of BPj only. The 
magnetostatic energy of such a Bloch point has the form 

1 

g™ = i J r 2 [7 + 4 cos 7 (r) + 4 cos 2 7 (r)] dr. (17b) 



From Eq. ( lc ) we obtain that the Bloch point interaction 



with magnetic field can be expressed as follows 



1 

S l = -26 J r 3 cos7(r)dr. 



(17c) 



By minimizing the total energy, Sg/Sj — 0, we obtain 
that the equilibrium distribution 7(r) is a solution of the 
following nonlinear differential equation 

d 2 7 2ed 7 1 2 

e-r-5- H ; — I- 7 sin 7 + - sin 27 — or sin 7 = (18) 

dr^ r dr 5 5 

augmented by boundary conditions of the form 

= 0. (19) 
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FIG. 3: (Color online.) Reduced rotation angle g £ (r), 
see ( |20| for different field intensities and e — 0.05: 
analytical result ( |Al[) (s olid curve) and numerical 
solution of (|18|) (dashed curves). 



In the case of weak fields one can linearize Eq. ( 18 1 in 



the vicinity of spatially uniform solution (13) and obtain 
that 



70") « 70 + bg e (r), \b\ < 1. 



(20) 



An explicit form of the function g £ (r) is calculated in 
Appendix [A] The comparison with numerical solution of 
Eq. ( |I8[ ) shows a quite good agreement up to relatively 
strong fields (b < 1), see Fig. [3] 

Another limiting case is realized in the case of strong 
magnetic fields when the Bloch point magnetization is 
parallel to the external field. In this case the rotation 
angle is 7 = (mod n). 

To describe the behavior of the Bloch point in a criti- 
cal region b ~ b c where the spatially non-uniform distri- 
bution transforms to the spatially uniform one, we use a 
variational approach with a two-harmonics trial function 
7(7") « ao + ai cos7rr. Near the critical point ao, ai <C 1. 
We expand the total energy in a Taylor series up to the 
fourth order with respect to ao and to the second order 
with respect to a%. By excluding a\ and keeping terms 
not higher than a^, we get 



^(7) « £'o+p2(b,e)al+p4(b,s)ao, 



(21) 



The energy (21) as a function of ao has a double-well 
shape (p2 (b, £j< 0) for b < b c with the critical magnetic 
field b c given by 

b c (e) w 1.8 - 21.6e + ^0.4 - 20.2e + 467e 2 (22) 

In the critical region when < b c (e) — b <C b c (e) , 

a (b)^a(e)^b c (e)-b (23) 



For b > 6 C , P2 > and the function Eq. (21) has a 
minimum for ao = 0. It corresponds to 7 = 0. Numer- 
ical integration of Eq. (18) for e — 0.05 shows that the 



phase transition occurs when b c « 1.473, see Fig. [4j It 
agrees well with the value 6 C (0.05) ~ 1.465 obtained from 
Eq. (p2l. The critical behavior predicted by |23l is also 



confirmed by our numerical simulations (see Fig. 4a) 



III. 



NUMERICAL STUDY OF THE BLOCH 
POINT STRUCTURE 



In order to check analytical results about Bloch 
point structure, we performed simulations using in-house 
developed spin-lattice simulator SLaSiP^l that solves 
Landau-Lifshitz-Gilbert equation in terms of spins 



dSn 
dt 



dJt? 



dS n 



dS n 
dt 



where ffl is a lattice Hamiltonian of the classical ferro- 
magnet: 



(n,5) 



2/4 £ 



(S n ■ gfe) _ 3 (S„ • r nk )(S k ■ r nk ) 



r n k- 



(24) 



Here S n is a classical spin vector with fixed length S 
in units of action on the site n of a three-dimensional 
cubic lattice with lattice constant ao, J is the exchange 
integral, /ib is Bohr magneton, r nk is the radius-vector 
between n-th and fc-th nodes, 77 is a damping param- 
eter, H is external magnetic field and 5 runs over six 
nearest neighbors. Integration is performed by modified 
4-5 order Runge-Kutta-Fehlberg method (RKF45) and 
free spins on the surface of the sample.^ 

Numerically we checked the Bloch point structure, 
given by the radial-dependent Ansatz (fl6|). by model- 



ing spherically-shaped sample with diameter 2R = 35ao 
(such a sample consists of 24464 nodes with nonzero 
spin), and exchange length £ = 3.95a (e = 0.05). In 
order to stabilize the Bloch point we applied the gra- 
dient magnetic field with b = 1.0. By modeling the 
overdamped dynamics we observed that the Bloch point 
structure quickly relaxes to the state similar to one, 
given by (16): The polar Bloch point angle 9(r) does 
not deviate from i9 within the accuracy 0.099. The az- 



imuthal angle is well also well-described by ( 16 ) with the 
radial-dependent rotation angle 7(7"), see Fig. |5| Simula- 
tions were performed for crystallographic directions [111] 
{■d = tt/4) and [110] (1? w tt/2, the plane is shifted by 
z = — 0.5ao from the origin). One can see from Fig.[5]that 
numerical data are well confirmed by analytical curve 
7(7-), calculated as numerical solution of (18). 

To validate our theory we performed also direct stabil- 
ity check. Numerically we check the stability of the Bloch 
point against the shift of its position. We start simula- 
tions with the Bloch point state using Ansatz-function 
(16), which is shifted along £-axis by Az = —2ao. We 
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FIG. 4: (Color online.) Bloch point under the action of the gradient field, (a): Rotation angle vs field intensity b 
near the critical field b c ps 1.473 from numerical solution of Eq. (18) (blue curves) and theoretical estimation by 
Eq. (231 (red curves) with e = 0.05. Solid lines correspond to the rotation angle 7(0) and dashed line to 7(1). (b): 
Phase diagram for solutions of Eq. (18). The upper (hedgehog) phase correspond to the solution 7 = 0, the lower 



(nonhomogeneous) one to the radial-dependent Bloch point with 7(7*). Dashed lines correspond to analytical result 

for the critical field b c = 1.465 for e = 0.05, see text. 
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FIG. 5: (Color online.) Radial dependence of rotation 
angle 7 in spherical particle. Line: numerical 
integration of ( |18[ ). Symbols: SLaSi simulations for 
crystallographic directions [110] and [111]. Parameters 
are the same as in Fig. [2] 



also apply j(r,t = 0) = 3° in order to break the sym- 
metry. For rapid relaxation we used in most of simu- 
lations the overdamped regime (the damping parame- 
ter T) — 0.5). We checked the shift of the Bloch point 
by controlling the total spin projections: only for the 
Bloch point, situated at the sample origin, the total spin 

ctot ctot Ctot pi 

°x — °y — °z — u - 

The temporal evolution of initially shifted Bloch point 
is presented in Fig. [6] for the Bloch point sample with 
2R = 35ao (24 456 nodes) in applied field with 6=1, 



see also the supplementary videoSS. Originally the Bloch 
point was shifted down from the origin which corresponds 
to 5* ot > 0, see inset (a). During the evolution a number 
of magnons are generated, inset (b). After quick damp- 
ing of oscillations, the micromagnetic singularity goes to 
the sample origin, see inset (c). The relaxation process 
consists of two parts: (i) The rotation angle 7(7-) changes 
its value from initial uniform one to the final nonhomo- 



geneous state during a time r 7 



500w n 



(ii) The re- 



laxation of 5* ot component of total spin of the sample 
tooks approximately the same time. During all simula- 
tions time |S^ ot | « |^ ot | < 10" 11 . 



IV. CONCLUSION 

To summarize, we study the magnetization structure 
of the Bloch point. In spite of the fact that the Bloch 
point as a simplest 3D topological singularity was stud- 
ied durin g a long time, from the pioneer papers by 
Feldtkeller^ and Doring^, see also for review Refs. |3] 
and [lOl the probl em of t he Bloch point structure still 
causes discussions ! 16 * 30 * 3 -^ The point is that the most 
strong exchange interaction determines only the relative 
magnetization distribution accurate within the rotation 
angle 7. This rotation angle, which is determined by the 
magnetostatic interaction, is most questionable: its value 
is equal to 120° according to FeldtkellerEl, to 112.3° fol- 
lowing Doring^ and 113° following Elias and Verga^. 
We analyze the origin of all these results and calculated 
the equilibrium value, about 105°, see (12), which mini- 



mizes the total magnetostatic energy, not only the part 
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FIG. 6: (Color online.) Dynamics of total spin along z-axis of the sample. The Bloch point is initially shifted by 
Az = — 2ao from center of the sample. Insets show magnetization distribution in z = — 0.5ao and y — — 0.5ao planes 
in different times. Color bar indicates S z ,n for different lattice nodes. Applied field amplitude 6=1, other 

parameters are the same as in Fig. [2| 



of it. 

The next problem appears in modeling of the Bloch 
point. As is was discussed by Thiaville et al.^, the mod- 
eling of singularity is mesh-dependent. In particular, a 
mesh-friction effect and a strong mesh dependence of the 
switching field during the Bloch-point-mediated vortex 
switching process was detected using OOMMF micromag- 
netic simulations.^ The reason is that micromagnetic 
simulators consider the numerically discretized Landau- 
Lifstitz equation, which are valid in continuum theory. 
Since the Bloch point appears as a singularity of con- 
tinuum theory, it is always located between mesh points, 
and causes the mesh-dependent effects and therefore may 
be insufficient for describing near-field Bloch point dis- 
tribution. In contrast to this, spin-lattice simulations are 
free from these shortage. From the beginning we consider 
discrete spins, located on the cubic lattice, and their dy- 
namics is governed by the discrete versions of Landau- 
Lifshitz equations. The lattice Hamiltonian allows us to 
calculate the discrete energy of the Bloch point similar 
to the atomiclike calculations by Reinhardt^. 

Using in-house developed spin-lattice SLaSi 32 simula- 
tor we modeled the Bloch point state nanosphere and 
checked our analytical predictions about Bloch point 
structure. We stabilized the singularity inside the spher- 
ical particle by applied gradient magnetic field. The field 
causes the new type of Bloch point with radial-dependent 
rotation angle j(r). 
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Appendix A: Bloch point structure in a weak field 

We consider here the magnetization structure of a 
Bloch point under the action of weak magnetic field. One 
has to linearize Eq. ( 18 ) on the background of the unper- 



turbed rotation angle 70, see (20 1, which can be presented 
as follows: 



2v5e 1 /3 

j(r) ~~/o + bg s (r), g e {r) = — — /(Ar), X= 2\~- 

Here the function /(£) satisfies the linearized version of 
Eq. (pi): 



dV 2d/ 
which can be easily integrated: 



f = t, 



/(o = c A si ^ + 2 cosh i^-e 



A 2 - 2A sinh A + 2 cosh A — 2 
A cosh A — sinh A 



(Al) 
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The graphics of the g £ (r) for e = 0.05 is presented in 
Fig. [3] together with numerical solution of Eq. ( 18 ) by 
shooting method. In spite of limitation of our analysis 
by the case of weak field, |6| <§C 1, the function g s (r) pro- 
vides a good approximation for the solution of nonlinear 
Eq. ( 18 ) up to to very strong fields b < 1 with a relative 
error \ ff(r) 

num - 7( r )thcor] /l(r) mini 

< 0.04. 
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The rotation angle in the Bloch point is essentially 
influenced by the exchange parameter e, see Fig. [7] In the 
limit case of small particle (s 3> 1) the role of exchange 
is dominant, which results in the constant angle 1700 = 
— \/l5/4 « —0.97. In the opposite case e -C 1, the role 
of magnetostatic interaction is enhanced and this leads 
to a nonhomogeneous rotational angle distribution. In 
the limiting case go(0) = and <?o(l) = —y/5/3 ~ —1.3. 
Such a limit case is realized in typical soft nanomagnets 
sized in some tens of nanometers. 



FIG. 7: (Color online.) Reduced rotation angle g e vs 
reduced exchange length e: at r = (solid curve) and 
r = 1 (dashed curve). 
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